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ABSTRACT 


A thermal model was developed to predict the thermal response of coated and uncoated 
tubes tested in a H 2 /0 2 rocket engine. Temperatures were predicted for traditional APS Zr0 2 -Y 2 0 3 
thermal barrier coatings, as well as APS and LPPS Zr 0 2 -Y 2 03 /NiCrAIY cermet coatings. Good 
agreement was observed between predicted and measured metal temperatures at locations near 
the tube surface or at the inner tube wall. 

The thermal model was also used to quantitatively examine the effect of various coating 
system parameters on the temperatures in the substrate and coating. Accordingly, the effect of 
the presence a metallic bond coat and the effect of radiation from the surface of the ceramic layer 
were examined. In addition, the effect of a variation in the values of the thermal conductivity of the 
ceramic layer was also investigated. It was shown that a variation in the thermal conductivity of 
the ceramic layer, on the order of that reported in the literature for plasma sprayed Zr0 2 -Y 2 0 3 
coatings, can result in temperature differences in the substrate greater than 100°C, a much 
greater effect than that due to the presence of a bond coat or radiation from the ceramic layer. 
The thermal model was also used to predict the thermal response of a coated rod in order to 
quantify the difference in the metal temperatures between the two substrate geometries in order 
to explain the previously-observed increased life of coatings on rods over that on tubes. It was 
shown that for the short duration testing in the rocket engine, the temperature in a tube could 
exceed that in a rod by more than 100°C. Lastly, a two-dimensional model was developed to 
evaluate the effect of tangential heat transfer around the tube and its impact on reducing the 
stagnation point temperature. It was also shown that tangential heat transfer does not 
significantly reduce the stagnation point temperature, thus allowing application of a simpler, one- 
dimensional model for comparing measured and predicted stagnation point temperatures. 
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Thermal Modelling of Various Thermal Barrier Coatings 
In A High Heat Flux Rocket Engine 


INTRODUCTION 

Previous testing (Ref. 1) of various thermal barrier coatings (TBC's) in a high-heat flux 
hydrogen/oxygen (H 2 /0 2 ) rocket engine at the NASA Lewis Research Center (NASA-LeRC) 
showed that these coatings exhibit significant potential to reduce metal temperatures during short 
duration testing. In addition, durability testing of numerous TBC's from a variety of sources 
showed that certain TBC's possess considerable durability in a high heat flux environment, 
although repeatability in the life of a coating was an issue. This durability study, carried out with 
both rod and tube substrates, found that, in general, coatings were more durable on rod 
substrates than on tubes. The increased durability was attributed to the greater thermal mass of 
the rod resulting in lower coating and substrate temperatures, although temperature 
measurements on rod substrates were not performed. 

The purpose of the present work was to model the thermal response of coated and 
uncoated tubes and rods tested in the rocket engine at NASA-LeRC. The coatings included 
traditional TBC's with a ceramic top coat and metallic bond coat, as well as two cermet coatings 
consisting of different amounts of ceramic and metal in the outer layer above a metallic bond coat. 
One cermet coating was applied by low pressure plasma spraying (LPPS) while the second was 
applied by air plasma spraying (APS). A thermal model was developed to predict temperatures 
for both tube and rod geometries in both the coated and uncoated conditions. These predictions 
were compared with metal temperatures measured just beneath the outer surface or at the inner 
wall of the tube (Ref. 1 ). Use was made of previously-measured gas temperatures and heat 
transfer coefficients for tubes located at the throat of the rocket engine (Ref. 2). 

The thermal model was also used to examine the effect of the variation in reported values 
of the thermal conductivity of plasma sprayed Zr0 2 -Y 2 0 3 ceramic coatings, the effect of the 
presence of a metallic bond coat, the effect of tangential heat transfer in reducing the temperature 
at the stagnation point of the tube, and the effect of radiation from the surface of the ceramic 
coating. The predicted effectiveness of the various coatings in reducing the thermal shock to 
SSME-HPFTP blades has been reported elsewhere (Ref. 3) 

Background to Previous Experimental Results 

Test Environment . Test tubes or rods were positioned at the throat of the gaseous H 2 /O z rocket 
engine. The samples were held in a copper holder such that all gas flow was perpendicular to the 
cylindrical axis of the samples. High pressures in the combustion chamber resulted in Mach 1 
exhaust gas velocities slightly upstream from the throat plane (the plane passing through the 
cylindrical axis of each of the samples). Five tubes or rods were held in the copper holder, 
however, gas temperature measurements indicated significant temperature differences between 
the five test positions (Ref. 2). Hence, ail durability testing of the coated tubes and rods had been 
performed in the center position. Metal temperatures in the tubes were only measured at the 
midpoint (lengthwise) of the tube where the gas temperature had been characterized in greatest 
detail (Ref 2). Gas temperatures at this position were relatively constant around the tube from the 
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stagnation point (0=0°, facing the injector) through the throat plane (0=90°), but decreased on the 
exhaust side of the tube opposite the injector (0=180°). Typical gas temperatures for the 
stagnation point and on the exhaust side of the tube are shown in Figure 1. These relatively low 
gas temperatures were achieved by maintaining a hydrogen-rich environment in the rocket 
chamber with an 0 2 to H 2 mass ratio of approximately 1 .2 to 1 .4 (stoichiometric combustion at a 
mass ratio near 8 produces gas temperatures in excess of 3000°C). The decrease in gas 
temperature on the exhaust side of the tube is likely the result of entrained ambient air due to a 
flow separation which occurs near 0=±13O° (Ref 2). Maximum heat transfer coefficients h c(max) 
around the tube (Ref 2) are shown in Fig 3a. The decrease in h^^ on the exhaust side of the 
tube is again due to the flow separation. The heat transfer coefficients also varied with time 
during engine operation due to changing pressures and purge gases in the rocket chamber (Ref 
2). The coefficients were at a maximum for 1 .2 seconds during the high pressure, high 
temperature portion of the test cycle (5-6.2 seconds, Fig 1), were lower for 0.4 seconds during a 
brief H 2 purge when the oxygen valve was closed and the pressure decreased, and were lower 
still for a final N 2 purge at approximately atmospheric pressure. The time dependence of h c was 
taken into account with a scaling factor k, shown in Fig 3b. Operation of the rocket engine and 
details of the measurements of the gas temperature and metal temperatures in the tubes are 
discussed in detail elsewhere (Refs. 2,4). 

Substrates and Coatings . The substrate material for most tubes and rods was directionally 
solidified (DS) Mar-M246+Hf 1 . One tube of conventionally cast B-1900 superalloy 2 was also 
coated and tested. Both tubes and rods had an outer diameter of 0.953 cm (0.375 in) and the 
tubes had an inner diameter of 0.648 cm (0.255 in). Temperatures were measured with 
thermocouples located just beneath the outer tube surface or at the inner tube wall. The 
instrumented tubes were positioned so that the thermocouples were nearest the injector (0=0, see 
Fig 1). Details of the thermocoupling procedures and techniques are given elsewhere (Ref 2). 
The tubes were coated with either traditional plasma-sprayed Zr0 2 -Y 2 0 3 TBC’s or with cermet 
coatings. The traditional TBC’s consisted of an APS Ni-Cr-AI-Y bond coat below a Zr0 2 -7Y 2 0 3 
top coat. Two thicknesses of coating were tested in either the as-sprayed condition or after a 
smoothing treatment with SiC paper. Two different cermet coatings were also tested. The bond 
coat and cermet layer of one coating was applied by LPPS, while both layers of the second 
coating were applied by APS. The LPPS cermet contained approximately 47% Zr0 2 -8Y 2 0 3 in a 
NiCrAlY matrix whereas the APS cermet contained only 33% of the same ceramic. Coating 
designations and thicknesses for both the traditional TBC’s and the cermet coatings are given in 
Table I. Densities of the different coating materials are given in Table II. Further details of the 
coatings, including compositions, surface roughness, plasma spraying parameters for the 
traditional TBC’s, and cross-sectional microstructures, are given elsewhere (Ref. 4) 


THERMAL MODELLING ANALYSIS AND RESULTS 
Thermal Model Development and Analysis . Two finite-difference thermal models were developed 


1 The nominal composition of DS Mar-M246+Hf is 

Ni- 1 0Co-9Cr- 1 0W-5.5Al-2.5Mo- 1 ,5Ta-l ,5Ti-0. 1 5C-0.0 1 5B-0.05Zr-l .75Hf weight percent. 

1 The nominal composition of B-1900 superalloy isNi-10Co-8Cr-6Mo-6Al-4Ta-lTi-0.1C-0.1Zr-0.015B weight 
percent 
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to simulate heat transport in coated and uncoated tubes and rods during testing in the rocket 
engine. The first thermal model considered convective heat transfer over the entire outer surface 
of the tube from the stagnation point (0=0°) to the exhaust side of the tube (0=1 80°) and 
accounted for both radial and tangential conductive heat transfer within the tube wall (2-D model). 
The second thermal model simulated convective heat transfer at the stagnation point and 
accounted for conductive heat transfer in a radial direction only (1-D model). In both models, the 
appropriate heat transfer equations for conduction, convection and radiation were solved by the 
numerical computer program known as SINDA (Systems Improved Numerical Differencing 
Analyzer) (Ref. 5). The SINDA program requires, as a minimum, the thermal conductivity (K), the 
specific heat (Cp), and the density (p) for all materials transferring heat by conduction. The 
program also requires a grid to be defined for the geometry of interest. If desired, the SINDA 
program allows time and temperature-dependent thermophysical properties as well as time- 
dependent gas temperatures and heat transfer coefficients (h c ). The program also accounts for 
heat loss from the surface by radiation. 

2-D Model 

The maximum effect of the tangential heat flow was evaluated with the 2-D model by 
simulating heat transfer to, and conduction around, an uncoated tube. This model consisted of 
110 nodes and equidistant grid spacing, as shown schematically in Figure 2. Only half of the 
cylinder cross section was modelled due to the symmetry in gas flow around the upper and lower 
halves of the cylinder. Based on gas temperature measurements made around the tubes (Ref. 

2), the higher gas temperature profile (0=0°, Fig 1) was used for the surface nodes between 0=0° 
to 0=150° and the lower temperature profile (0=180°, Fig 1) for the surface nodes between 0=150° 
to 0=180°. Values for h c at the 10 surface nodes were taken from Fig 3a. The time dependence 
of h c was incorporated by multiplying the value for hc( max) from Fig 3a with the time-dependent 
scaling factor, <|>, from Fig 3b. Hence, each of the surface nodes have the value for h c(max) shown 
in Fig 3a between the time of 5.2 to 6.2 seconds, approximately 25% of that value during the H z 
purge between 6.2 and 6.6 seconds, and approximately 5% of the maximum value during the N z 
purge, thereafter. The temperature dependence of Cp and K for the DS Mar-M246+Hf substrates 
(Ref. 6) is shown in Figures 4a and b, respectively, and the density p (Ref. 6) is given in Table II. 

Surface temperatures around the tube predicted with the 2-D thermal model after 0.3, 0.6 
and 1.2 seconds are shown in Figure 5a. As expected, the surface temperature at the stagnation 
point rises more rapidly than that on the exhaust side of the tube, rising to more than 80% of the 
maximum temperature after only 0.3 seconds. To quantitatively assess the effect of the tangential 
heat flow, the 2-D thermal model was also run with zero conductance between the nodes in the 
circumferential direction (i.e., radial heat flow only, no tangential flow). The time dependence of 
the predicted surface temperature at 0=0° and at 0=180°, with and without tangential heat flow, is 
shown in Figure 5b. There is less than a 10°C reduction in the stagnation point temperature 
(0=0°) when tangential heat transfer is taken into account. However, on the exhaust side of the 
tube (0=180°), heat flow around the tube to the exhaust side increases the temperature 
approximately 140°C. Hence, tangential heat flow does not significantly reduce the stagnation 
point temperature but does substantially increase the temperature on the exhaust side of the tube. 
Since all temperature measurements in the tubes were made at a location facing the injector (0=0, 
Fig 1), the simpler, 1-D thermal model was thereafter used to predict the stagnation point 
temperatures measured in the tubes. 

1-D Model 
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Uncoated Tube : The 1-D thermal model consisted of a 45° wedge that contained 61 nodes for 
the tube and 86 nodes for the rod. For the tube, the near-surface grid spacing was 12.5 pm 
(0.0005 in) and the spacing near the inner wall was 50 pm (0.002 in). The near-surface grid 
spacing for the rod was the same as that for the tube increasing to 125 pm (0.005 in) at the center 
of the rod. The 1-D model is shown schematically in Figure 6. Thermophysical properties for the 
DS Mar-M246+Hf substrate (Cp, K and p) were used as with the 2-D model (i.e., Figures 4a, b and 
Table II). The gas temperature (Fig 1) and heat transfer coefficient (Fig 3a) appropriate for the 
stagnation point (0=0°) were used to conduct heat to the surface node. The maximum value of the 
heat transfer coefficient (33 kW/m 2 K) was scaled with time according to the factor shown in Fig 
3b. The convergence of the numerical solution (i.e., an adequate number of nodes) was verified 
by decreasing the number of nodes from 61 to 16. Temperature profiles through the wall of the 
tube differed by less than 1°C when 30 nodes were used. However, when the number of nodes 
was reduced to 16, surface temperatures were up to 9° higher while inner wall temperatures were 
up to 41 °C lower. Hence, the grid spacing resulting from using 30 or more nodes was sufficient to 
produce a convergent solution. Predicted and measured metal temperatures for the uncoated 
Mar-M246+Hf tube are shown in Fig 7. The agreement is good considering the uncertainties of 
the gas temperature measurements (Ref 2). 

Coated Tube : Prior to predicting temperature profiles for the coated tubes tested in the rocket 
chamber, the 1-D model was used to examine various aspects associated with modelling the heat 
transfer to coated tubes. Since there is considerable variation in the reported values for the 
thermal conductivity of APS Zr0 2 -Y 2 0 3 ceramic coatings, the effect of this variation on predicted 
temperatures in the ceramic layer and the metal substrate was examined. In addition, the need to 
account for the presence of a metallic bond coat, and the need to account for radiation from the 
surface, were also investigated. Thermophysical properties for each of the coatings, whether 
metallic, cermet or ceramic, were assigned to the outer grids of the 1-D model at a maximum grid 
spacing of 12.5 pm (0.0005 in). Coating designations, densities and reference sources for the 
various coatings are given in Table II. 

Several reported values of the temperature dependence of the specific heat (Cp) and the 
thermal conductivity (K) for APS Zr0 2 -Y 2 0 3 are shown in Figures 8a and 8b. The effect of the 
variation in the thermal conductivity values was examined by predicting temperature profiles for a 
180 pm (0.007 in) thick ceramic layer on a Mar-M246+Hf tube. For the first case, (I), lower values 
for K were examined using the data from Ref 10 (Figure 8b). For the second case (II), the higher 
values for K from Ref 7 were used (Figure 8b). The nearly identical data from Refs 7 and 8 were 
used for Cp (Figure 8a). All thermophysical data for the Mar-M246+Hf substrate was as 
described above for the uncoated tube. The predicted temperatures for the surface of the ceramic 
and at the ceramic/metal interface for these two cases are shown in Figure 9. The variation in the 
values of the thermal conductivity amounts to less than a 50°C difference in temperature at the 
ceramic surface but a 180°C difference at the ceramic/metal interface. Hence, the range in the 
reported values for K for APS Zr0 2 -Y 2 0 3 coatings result in significant temperature differences at 
the ceramic/metal interface but only moderate differences at the surface of the ceramic. 

The TBC's modelled in this study contained a metallic bond coat beneath the cermet or 
ceramic layer. These bond coats were 25-75 pm (0.001-0.003 in) thick and were either plasma 
sprayed in air (APS) or at reduced oxygen pressures (LPPS). Values for Cp and K for APS and 
LPPS metallic coatings are shown in Figures 4a and b, respectively. Most of the Cp data are 
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closely grouped, however, the thermal conductivity of the APS coatings is generally much lower 
than that for the LPPS coatings. The lower values for the APS coating can generally be 
attributed to the amount and structure of the oxide present in these coatings. The maximum 
effect of a bond coat was evaluated by predicting substrate temperatures for a TBC with a 100 
p.m (0.004 in) APS bond coat below a 180 jam (0.007 in) ceramic top coat. Predictions were 
made with higher (Case I) and lower (Case II) thermal conductivity values for APS coatings. At 
the higher temperatures, the values for K differ by a factor of two (see data from Refs. 7 and 1 1 
in Fig 4b). Similar values for Cp for the bond coat and identical values for K and Cp for the 
ceramic layer were used. Figure 10 shows the predicted temperature just below the bond 
coat/substrate interface for the two cases as well as the predicted temperature at the same 
location without a bond coat. For Case I with the higher conductivity bond coat, the 
temperature in the substrate was reduced by only 10°C while for Case II with the lower 
conductivity bond coat, the temperature in the substrate was reduced by approximately twice 
this amount. Temperatures increased in the ceramic layer, but by smaller amounts than the 
decreases in the substrate. Since the thermal conductivity for LPPS coatings is greater than 
that for the APS coatings and similar to, or greater, than that for the substrate (Figure 4b), an 
LPPS bond coat would be expected to have a negligible effect on the temperatures in the 
substrate or the ceramic coating. Hence, the bond coat produces only a small effect on the 
predicted temperatures in either the ceramic layer or the substrate. Consequently, bond coats 
were not included in the thermal model when predicting temperatures for the coated tubes 
tested in the rocket chamber. 

Radiation from the ceramic coating reduces the temperature of the surface. Liebert 
(Ref. 13) has reported correlation equations for the emissivity for several thicknesses of APS 
Zr0 2 -Y 2 0 3 TBC's in the temperature range of 27°C to 2500°C. For a 180 jam (0.007 in) thick 
coating, the work by Liebert suggests a linearly decreasing value of the emissivity from 0.96 at 
27°C to 0.46 at 1500°C. These emissivity values were used in the thermal model to evaluate 
the reduction in surface temperature due to radiative heat losses. Accounting for heat loss by 
radiation amounted to less than a 10°C decrease at the surface of the ceramic coating and an 
almost negligible decrease at the ceramic/metal interface. Hence, radiative heat losses from 
the ceramic coating were also ignored when predicting temperatures for the coated tubes 
tested in the rocket chamber. 

Comparison of Measured and Predicted Metal Temperatures in TBC-Coated Tubes . The 1-D 
tube model was used to predict temperatures in the traditional TBC's and the APS and LPPS 
cermet coatings. Based on the previous modelling results, a metallic bond coat and radiative heat 
losses were not included in the model. The gas temperature and heat transfer coefficient were as 
previously used for the 1-D model. Thermophysical data for the DS Mar-M246+Hf substrate were 
also used as previously discussed. Median values for the thermal conductivity of the Zr0 2 -Y 2 0 3 
coatings were defined and used (Figure 8b). Values for Cp from Reference 8 (Figure 8a) were 
taken as representative of the Zr0 2 -Y 2 0 3 coatings and a value of 4.8 gm/cm 3 was taken as a 
typical density (Table II). For the B-1900 substrate, values for Cp, K and the density of Mar- 
M246+Hf were used because of the similarity in values between the two Ni-base alloys (Ref. 14). 
The measured and predicted temperatures just below the ceramic/metal interface for the thick 
ceramic coatings (200 and 250 |im thickness) are shown in Figure 11a and for the thinner 
coatings (100 and 140 p,m thickness) in Figure 11b. The measured and predicted inner wall 
temperatures for the thick ceramic coating (250 pm thickness) on the B-1900 substrate are shown 
in Figure 11c. The decreasing temperatures associated with the H 2 purge cooling period (6.2-6.5 
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seconds) is predicted somewhat sooner than actually observed. This difference in the onset of 
cooling may indicate the finite time required for the 0 2 valve to close and for the H 2 purge to 
actually begin. However, in considering the variation in the reported values for the thermal 
conductivity of the ceramic layer, the agreement between the measured and predicted 
temperatures is quite good. 

The temperature dependence of Cp and K for the cermet coatings is given in Figures 4a 
and b. Average values for the coating thicknesses of 64 jam (0.0025 in) for the LPPS and 180 
pm (0.007 in) for the APS cermet coatings (see Table I) were used. Measured and predicted 
inner wall temperatures for these two coatings are shown in Figure 12. Again, the agreement is 
quite reasonable considering the number of potential sources of uncertainty. 

Predicted Metal Temperatures for Traditional TBC-Coated Tubes vs. Rods . The difference in the 
thermal response between a tube and rod was examined by predicting temperature profiles for 
both substrate geometries with 100 pm and 200 pm (0.004 and 0.008 in) thick ceramic coatings. 
The gas temperature and heat transfer coefficient were those discussed above for the 1-D model. 
Thermophysical properties and densities for the Mar-M246+Hf substrate and for the ceramic layer 
were as used in the proceeding section. The temperatures predicted on the surface of the 
ceramic layer and at the ceramic/metal interface for both the tube and rod geometry are shown in 
Figure 13a. There is only a small difference in the surface temperature of the ceramic between 
the tube and rod where, after 1 .2 seconds, the rod was 34° and 1 6°C lower than the tube for the 
thinner and thicker coating, respectively. However, at the ceramic/metal interface, the 
temperature for the rod was 160°C and 145°C lower than the tube for the thinner and thicker 
coatings, respectively. These differences in temperature between the tube and rod reflect the 
greater thermal mass of the rod. Hence, for the given gas temperature profile, the metal surface 
temperature of the rod is approximately 1 50°C lower than that for the tube. 

Figure 13a also highlights the importance of coating thickness on the substrate 
temperature. There is little temperature difference at the ceramic surface between the 100 and 
200 pm coatings for either the tube or rod ( a maximum of 60°C after 0.3 sec). However, at the 
ceramic/metal interface, the temperature difference between the two coating thicknesses is 224°C 
(+8°C depending on the tube or rod geometry). The predicted surface temperature of an uncoated 
tube after 1 .2 seconds, is approximately 1370°C (Fig 7). Hence, Figure 13 shows that the 
predicted benefit of a 100 pm (0.004 in) thick traditional TBC coating for the test conditions in this 
study is slightly under 300°C and the benefit of a 200 pm (0.008 in) thick traditional TBC is better 
than 500°C. As expected, the predicted temperature benefit is not directly proportional to the TBC 
thickness. 

Radial temperature profiles through a 200 pm (0.008 in) thick traditional TBC coating and 
into a rod and tube substrate are shown in Figure 13b. Simply considering the area under the 
temperature profile for each substrate shows that the rod acts as a greater heat sink than the 
tube. It appears that the temperature across the ceramic coating decreases linearly. However, a 
close examination of the temperatures indicated a slight concave downward curvature consistent 
with a nearly-constant surface temperature (see Figure 13a) and a positive temperature 
dependence of the thermal conductivity for the ceramic (Figure 8b). As the surface temperature 
rises at short times (<0.1 seconds), the temperatures in the ceramic layer did exhibit a slightly 
concave upward curvature, as expected during the early heat-up period. 
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DISCUSSION 


Thermal modelling allowed investigation of various aspects associated with heat transfer 
to coated tubes and rods (e.g., tube vs. rod substrates, the effect of the presence of a bond coat 
and the effect of radiation). In addition, thermal modelling also allowed the effect of the variation 
in the reported values of the thermophysical properties (specifically the thermal conductivity) to be 
examined and quantified. For a 100 pm (0.004 in) APS metallic bond coat, where the reported 
values for K differ by more than a factor of two, it was shown that the substrate temperatures were 
reduced by only 10° to 20°C while temperatures in the ceramic increased even less. This limited 
effect of the bond coat is expected by realizing that the overall thermal resistance of the TBC 
consists of the sum of the thermal resistance for each layer (ceramic and bond coat). Since the 
thermal resistance is inversely proportional to the thermal conductivity, the thermal resistance of 
the ceramic far outweighs that for the bond coat for similar layer thicknesses. The relatively high 
thermal conductivity of the LPPS bond coats, which is similar to that of the Mar-M246+Hf 
substrate, results in very little thermal resistance and produces little effect on temperatures in the 
substrate. 

The variation (approximately ±50%) in the reported values for the thermal conductivity of 
the APS Zr0 2 -Y 2 0 3 ceramic layer has a significant effect on the predicted tube substrate 
temperatures, and to a lesser extent, the temperature in the ceramic layer itself. It is not known 
whether the differences in the values for K are solely due to variations in the ceramic material 
(e.g., extent of porosity and microcracking, etc.) or to differences in the measuring techniques. It 
is not likely that small differences in the composition of the APS Zr0 2 -Y 2 0 3 ceramic can account 
for these differences in K. However, the good agreement between the measured and predicted 
temperatures for the traditional ZrO z -Y 2 0 3 TBC's suggests that thermal modelling can be 
successfully performed using median values for K. 

Modelling the thermal response of a coated tube and rod allowed the temperature 
difference due to the greater thermal mass of the rod to be quantified. In past durability testing of 
TBC's in the rocket engine (Ref. 1), coated rods were tested at gas temperatures estimated to be 
1 50-200°C greater than that for the coated tubes. Even at this higher temperature, the coating 
lives on the rods were generally greater than those for the tubes. The results presented in Figure 
13 show that for the same test conditions, the temperature at the ceramic/metal interface is 145- 
160°C lower for the rod than for the tube. Temperature profiles in a coated rod were also 
predicted for higher gas temperatures. For a 100 pm (0.004 in) thick TBC, these predictions show 
that the gas temperature for the coated rod must be 250°C greater than that for the coated tube in 
order to produce similar metal surface temperatures for both substrate geometries. Likewise, for 
a 200 pm (0.008 in) thick coating, the gas temperature for the coated rod must be increased 
310°C over that for the coated tube to produce similar metal surface temperatures. 

Although lower substrate temperatures in the rod result in a higher thermal gradient in the 
ceramic layer (see inset in Fig 13b) and presumably a greater thermal gradient-induced stress 
within the ceramic, most TBC failures are attributed to thermally induced stresses produced by a 
mismatch in the thermal expansion coefficients (CTE) between the ceramic and the metallic bond 
coat/substrate. It is suggested in the following discussion that higher temperatures in the tubes 
eventually produce greater stresses across the ceramic/metal interface resulting in shorter TBC 
lives in the tubes. Generally, a "stress free" temperature is defined as the temperature of the 
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substrate during application of the ceramic layer. Cooling a coated component below this stress 
free temperature results in a compressive load on the ceramic and tensile load on the metallic 
bond coat. Testing at temperatures above the stress free temperature results in a tensile load on 
the ceramic layer while the bond coat, constrained by the ceramic layer, is compressively loaded. 
High temperatures at the ceramic/metal interface, which occur during testing, can allow the 
relatively soft metallic bond coat to plastically deform by thermally-activated processes to reduce 
the compressive stress. The consequence of this relaxation by the bond coat is that the stress 
free temperature increases. Higher bond coat temperatures allow the bond coat to deform more 
quickly and easily and result in a higher stress free temperature. Upon cooling, this increase in 
the stress free temperature causes a greater compressive stress on the ceramic layer, which can 
result in cracking and spalling of the layer. Hence, it is suggested that the shorter TBC lives on 
tubes in comparison to rods is directly related to the higher metal temperatures encountered by 
the tubes, given similar TBC thicknesses for the tube and rod. It should be noted that this 
argument assumes that compressive loading within the ceramic due to rapid heating is insufficient 
to cause ceramic spallation (i.e., spalling occurs only on cooldown and not on heatup), and also 
that sufficient time is spent above the stress free temperature to allow the thermally-activated 
processes to occur. 

The good agreement between the predicted and measured temperatures for each of the 
various TBC's verifies the ability of the thermal model to simulate heat transfer to coated throat 
tubes during short thermal excursions. Good agreement was observed at locations near the 
ceramic/metal interface and at the inner wall of the tube. The slight time shift (0. 1 to 0.2 seconds) 
between the measured and predicted maximum temperatures at the ceramic/metal interface 
(Figures 11a and b), likely reflects the finite time required for the oxygen valve to close and for the 
hot gas to be expelled from the combustion chamber. However, the good overall agreement 
between predicted and measured temperatures justifies the use of thermal models to predict heat 
transfer during short thermal excursions using median thermophysical property data for the 
various coatings. 


SUMMARY AND CONCLUSIONS 

A thermal model was used to successfully predict metal temperatures in coated and 
uncoated tubes tested in an H 2 /0 2 rocket engine. Traditional APS Zr0 2 -Y 2 0 3 TBC’s, as well as 
APS and LPPS cermet coatings, were modelled. Good agreement was observed between 
predicted and measured temperatures in the metal near the tube surface for the uncoated tube 
and at the inner wall or near the coating/metal interface for the coated tubes. 

The thermal model was also used to examine the effect of the variation in the reported 
values of the thermal conductivity of plasma sprayed Zr0 2 -Y 2 0 3 ceramic coatings, the effect of a 
100 nm (0.004 in) thick metallic bond coat, the effect of tangential heat transfer around the tube, 
and the effect of radiation from the surface of the ceramic coating. It was shown that for the short 
duration rocket engine testing modelled in this study, radiation from the ceramic coating, 
tangential heat transfer within the tube wall, or the presence of a 1 00 (0.004 in) bond coat do 

not significantly affect the predicted metal temperatures. However, the difference in the reported 
values for the thermal conductivity of traditional APS Zr0 2 -Y 2 0 3 TBC's can amount to significantly 
different temperatures (>100°C) predicted in the tube. The thermal model was also used to 
compare the thermal response of a coated tube and rod. The results indicated substrate 


9 



temperatures approximately 1 50°C lower in the rod than in the tube. The increased life of 
coatings on rods over tubes tested in the rocket engine was directly attributed to this temperature 
difference. 
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FIGURE CAPTIONS 


Figure 1 . Representative gas temperature profiles at the stagnation point (0=0°) and at the 
exhaust side (0=180°) of the test position (Ref. 2). 

Figure 2. Schematic view of the 2-D thermal model of a tube incorporating radial and tangential 
heat transfer. Arrows indicate points of convective heat transfer. 

Figure 3. (a) Angular dependence of the maximum heat transfer coefficient around the tube, and 
(b) time dependence of the scaling factor <(> used to scale the heat transfer coefficient (h c ) to 
account for changes in pressure and purge gasses (Ref. 2). Gas temperature profile is also 
shown for comparison. 

Figure 4. Temperature dependence of the (a) specific heat (Cp) and (b) thermal conductivity (K) 
for the APS and LPPS metallic coatings, the cermet coatings, and the DS Mar-M246+Hf 
substrate. 

Figure 5. Predicted temperature profiles for an uncoated tube using the 2-D model (a) angular 
temperature distribution on the surface of the tube at various times with tangential heat flow, (b) 
time dependence of the surface temperature at 0=0 (stagnation point) and at 0=180° (exhaust 
side), with and without tangential heat flow. 

Figure 6. Schematic view of the 1-D thermal model of the rod and tube for heat transfer to the 
stagnation point and radial heat transfer. Arrow indicates the point of convective heat transfer. 

Figure 7. Measured and predicted metal temperatures for an uncoated tube at a location just 
below the outer tube surface (0=0). 

Figure 8. Temperature dependence of the (a) specific heat (Cp), and (b) thermal conductivity (K) 
for the Zr0 2 -Y 2 0 3 ceramic coatings. Specific heat for DS Mar-M246+Hf shown for comparison in 
(a). 

Figure 9. Predicted temperature profiles showing the effect of the variation in the values of the 
thermal conductivity for the Zr0 2 -Y 2 0 3 ceramic coatings. Case I used the lower values for K 
(Ref. 10) and Case II used the higher values for K (Ref 7) shown in Figure 8b. 

Figure 10. Predicted temperature profiles showing the effect of a 100 pm (0.004 in) APS metallic 
bond coat below a ZrO z -' Y 2 0 3 ceramic top coat. Predicted temperatures are for a location just 
below the bond coat/substrate interface. Case I used the higher K values (Ref. 11) and Case II 
used the lower K values (Ref 7) shown in Figure 4b. 

Figure 1 1 . Measured and predicted temperatures for the traditional Zr0 2 -Y 2 0 3 TBC's (a) near- 
surface temperatures for TK1 and TK2 with the thick ceramic coatings, (b) near-surface 
temperatures for TNI and TN2 with the thin ceramic coatings, and (c) inner wall temperatures for 
TK3, also with the thick ceramic coating. 
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Figure 12. Measured and predicted inner wall temperatures for the (a) LPPS and (b) APS cermet 
coated tubes. 

Figure 13. Predicted temperature profiles for a 100 and 200 jim (0.004 and 0.008 in) traditional 
Zr 02 -Y 2 0 3 TBC on a tube and rod substrate, (a) Time dependence of the temperature at the 
ceramic surface and at the ceramic/metal interface, and (b) radial temperature profiles in the 
ceramic and metallic substrate for the 200 (j.m (0.008 in) thick TBC's. 
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TABLE I 

Coating Designations, Thickness and Oxide in Top Coat 


Coating Designation 

Bond Coat 
pm (0.001 in) 

Top Coat 
pm (0.001 in) 

Percent Oxide 
in Top Coat 

Traditional TBC’s 3 




TNI 

25-50(1-2) 

100-125 (4-5) 

100 

TN2 

unmeasured prior to sanding (see below) 

100 

TK1 

25-50 (1-2) 

200-225 (7.9-8.8) 

100 

TK2 

unmeasured prior to sanding (see below) 

100 

TK3 d 

50-75 (2-3) 

225-275 (8.8-10.8) 

100 

After sanding 




TN2 

25-75(1-3) 

100-140(4-5.5) 

100 

TK2 

25-50 (1-2) 

200-250 (7.9-9.8) 

100 

Cermet Coating 0 




LPPS 

50 (2) 

50-75 (2-3) 

47 ±12 

APS 

40-60 (1. 6-2.4) 

150-200 (6-8) 

33 ± 10 


a APS NiCrAlY bond coat below APS Zr0 2 -7Y 2 0 3 ceramic top coat. TN refers 
to thinner ceramic top coat, TK refers to thicker ceramic top coat. 
b Measured from cross sections after testing 
c NiCrAlY bond coat below ZrO z -8Y 2 0 3 top coat 
d B-1900 superalloy substrate 
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TABLE II 

Designation and Density of Coatings and Substrate 


Designation 

Density 

(gm/cm 3 ) 

Reference 

Substrate 



DS Mar-M246+Hf 

8.6 

6 

APS Ceramic Coatings 



Zr0 2 -8Y 2 0 3 

3.47 

7 

Zr0 2 -(7-8)Y 2 0 3 

5.39 a 

8 


4.5 b 

9 

APS Metallic Coatings 



NiCrAI Y 

4.48 

7 

NiCrAI 

6.9 

8 

NiCr 

6.8 

11 

LPPS Metallic Coatings 



NiCrAlY 

5.48 

7 

NiCrAI 

7.48 c 

8 

Cermet Coatings 



LPPS 50/50CM 

4.55 

7 

APS 30/70CM 

4.01 

7 


a Representative values for several reported APS Zr0 2 -(7-8)Y 2 0 3 coatings 
b Based on reported 20% porosity 

c Average for three different LPPS NiCrAI based metallic coatings 
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Figure 1 . — Representative gas temperature profiles at 
the stagnation point (0 = 0°) and at the exhaust side 
(0 = 180°) of the test postion (Ref. 2). 
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Figure 2. — Schematic view of the 2-D thermal model 
of a tube incorporating radial and tangential heat 
transfer. Arrows indicate points of convective heat 
transfer. 
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Figure 3.— (a) Angular dependence of the maximum heat transfer coefficient around the tube, and (b) time depend- 
ence of the scaling factor <j> used to scale the heat transfer coefficient (h c ) to account for changes in pressure and 
purge gasses (Ref. 2). Gas temperature profile is also shown for comparison, 
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Figure 4. — Temperature dependence of the (a) specific heat (Cp) and (b) thermal conductivity (IQ for the APS and 
LPPS metallic coatings, the cermet coatings, and the DS Mar-M246+Hf substrate. 
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Figure 5. — Predicted temperature profiles for an uncoated tube using the 2-D model (a) angular temperature dis- 
tribution on the surface of the tube at various times with tangential heat flow, (b) time dependence of the surface 
temperature at 0 = 0° (stagnation point) and at 0 = 1 80° (exhaust side), with and without tangential heat flow. 
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Figure 6. — Schematic view of the 1 -D thermal model 
of the rod and tube for heat transfer to the stagnation 
point and radial heat transfer. Arrow indicates the 
point of convective heat transfer. 
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Figure 7. — Measured and predicted metal tempera- 
tures for an uncoated tube at a location just below 
the outer tube surface (0 = 0°). 
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Figure 8.— Temperature dependence of the (a) specific heat (Cp), and (b) thermal conductivity (K) for the Z1-O2-Y2O3 
ceramic coatings. 
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Figure 9.— Predicted temperature profiles showing the 
effect of the variation in the values of the thermal 
conductivity for the Z1O2-Y2O3 ceramic coatings. 
Case I used the lower values for K (Ref. 10) and 
Case II used the higher values for K ( Ref. 7) shown 
in Figure 8b. 
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Figure 10. — Predicted temperature profiles showing 
the effect of a 1 00 nm (0.004 in) APS metallic bond 
coat below a ZrO^^O^ ceramic top coat. Predict- 
ed temperatures are for a location just below the 
bond coat/substrate interface. Case I used the high' 
er K values (Ref. 1 1) and Case II used the lower K 
values (Ref. 7) shown in Figure 4b. 
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Figure 1 1 . — Measured and predicted temperatures for the traditional Z 1 O 2 -Y 2 O 3 TBC's (a) near-surface 
temperatures for TK1 and TK2 with the thick ceramic coatings, (b) near-surface temperatures for TNI 
and TN2 with the thin ceramic coatings, and (c) inner wall temperatures for TK3, also with the thick 
ceramic coating. 
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Figure 12. — Measured and predicted inner wall temper- 
atures for the LPPS and APS cermet coated tubes. 
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Figure 13. — Predicted temperature profiles for a 100 and 200 pm (0.004 and 0.008 in.) traditional Zr02~Y203 TBC on 
a tube and rod substrate, (a) Time dependence of the temperature at the ceramic surface and at the ceramic/metal 
interface, (b) Radial temperature profiles in the ceramic and metallic substrate for the 200 pm (0.008 in.) thick TBC's. 
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